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The kinetic roughening of a driven interface between three dimensional spin-up and spin-down 
domains in a model with non-conserved scalar order parameter and quenched disorder is studied 
O ' numerically within a discrete time dynamics at zero temperature. The exponents characterizing the 

\ morphology of the interface are obtained close to the depinning transition. 
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I. INTRODUCTION 



A variety of interface roughening models with quenched disorder have been studied recently (for an overview see 
jjj and references therein). These models have in common the occurrence of a so called depinning transition from a 
phase where the interface is trapped to a moving phase. To understand this transition and the morphology of the 
, — ' moving interface in the presence of quenched noise one may start with the continuum equation 

^ = t ,V 2 h + F + r ? (x,h(x,<)) (1) 



where h(x, t) denotes the interface profile, F a homogenous driving force and ??(x, h(x, t)) a noise term which is time 
, independent and depends only on the local position of the interface. The first term on the right hand side of this 
£h ' equation which is also known as Edwards- Wilkinson (EW) equation models a surface tension having a smoothing 
effect on the interface while the noise term roughens the interface. As has been shown by Bruinsma and Aeppli 
[0], Eq. (jl]) arises naturally in magnetic systems where the interface is between a domain with spins pointing in one 
direction and another domain with spins pointing in the opposite direction. Additionally, it has also been used for 
y—i , describing fluid flow in porous media Wa. 

Eq. ([!]) was modified by Kardar et al. || who added the nonlinear term -|(Vh) 2 with A proportional to the interface 
velocity. Such a term and many others arise when one tries to derive a continuum equation for the position of the 
interface from a Ginzburg-Landau-type free energy functional || . This approach can be traced back to Ref. JtJ where 
\ it was shown that the component of the local velocity normal to the domain wall is proportional to its local curvature. 
Expressing the position of the wall by a function h(x,i) assuming that there are no overhangs it is only a simple 
\Q mathematical exercise to derive from this the following equation of motion valid on length scales large compared to 
ON . the intrinsic width of the interface, 



c3 ■ 1 dh 



■3/2 



V 2 h + F + 77(x,h(x,i)), (2) 



where g = 1 + (Vh) 2 . The EW and the KPZ equation, respectively, follows by expanding in (Vh) 2 for nearly flat 
interfaces. 

There exists by now an extensive literature on these nonlinear equations as well as on various surface or interface 
models. Numerical studies of these models gave values for the characteristic critical exponents with rather large 
scattering. The question whether these numerical models fall into one of the two universality classes defined by the 
EW equation or the KPZ equation is therefore not always clear. 

In the present paper we add to this ongoing discussion a numerical study on the morphology of a domain wall in 
a three dimensional magnetic medium with quenched random fields which is driven by an external magnetic field. 
The magnetic medium is described by a Ginzburg-Landau-type energy functional and a zero temperature Langevin 
dynamics is imposed. We start from this semi-microscopic model since it is certainly more realistic than an equation 
of motion for the interface which is only obtained after quite a few more or less plausible approximations. 

It has been argued || that a flat domain wall just above the depinning transition belongs to the universality class of 
Eq. (||). As mentioned above to come to this conclusion a number of assumptions have to be made. The velocity of the 
domain wall must go to zero so that a KPZ-like quadratic term of kinetic origin can be neglected, multiplicative noise 
must be assumed to be irrelevant and the driving field must be replaced by an effective field which is the difference 
to the field strength at the depinning transition. Since these assumptions are of heuristic nature only it is of interest 
to study the underlying semi-microscopic model. 
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II. MODEL 



We study a model with non-conserved scalar order parameter (model A in the classification of Hohcnberg and 
Halperin |9|) with Langevin dynamics at zero temperature, 

9<f>(r,t) on 

7 dt d(t>(r,t)' [> 

with a relaxation time proportional to 7. Here thermal noise is neglected since it is believed to be irrelevant |n| at 
low temperatures. The Ginzburg-Landau type Hamiltonian TL is given by 

n = j dr(-^(r,t) 2 + ^(r ! t) 4 + ^(V^(r ! t)) 2 -(if + J B(r))^(r, i )) (4) 

where <j) denotes a scalar order parameter, H denotes a homogeneous driving magnetic field and B(r) a quenched 
random field. It is assumed that the random fields have zero mean and are uncorrelated in space and time. This 
model and the following discretization are discussed in more detail for the two dimensional or d = 1 + 1 case in a 
previous paper [TT| . 

The discretization of Eqs. (||) and (Q) for the computational simulation in the present case d = 2 + 1 is straight 
forward. It results in a set of difference equations, 



Si(t + At) = Si(r) + At 



(S 2 (t) - 1) Si(t) + ft'M - 65 i( r ) + h + h 



(5) 



which have to be iterated. The local magnetizations Si may be termed soft Ising spins at lattice point 1 = (x,y, z) 
with — 00 < Si < 00. The summation in Eq. (4) is over the 6 nearest neighbors S[ of Si. The driving field h and the 
quenched random fields bi are measured in units of J/6 2 , time r is measured in units of jS 2 / J and u = a5 2 / J = bS 2 / J 
whereas 6 denotes the lattice constant of the cubic lattice. Eq. (3) is iterated starting from a vertical flat initial 
interface where all spins on the left hand side of the interface located at x = are set to Si — 1 and all spins at the 
right hand side are set to S\ = — 1. In y- and z-direction periodic boundary conditions are assumed. The random 
fields bi are drawn with equal probability from an interval between —p and p with p = 0.6, u is chosen as u = 0.9 and 
the time constant as At = 0.1. 



III. RESULTS 



The interface location h(ro, t) is defined as that point at which the magnetization Si{t) as a function of x for fixed 
r o = {Vi z ) changes sign. For not too large values of p and h the function h(ro,T) is single valued since there are 
practically no overhangs or droplets. Of interest are the height correlation function 

C(r, t) = ([h(r +r,T)-h(r ,T)] 2 ) . (6) 

and the width of the interface, which is often called the roughness, 

w(L, t) = mv^r) - {Hr^r))) 2 ) 1 ' 2 (7) 

where the angular brackets denote an averaging over all interface sites and the over-bar the average over different 
realizations of the quenched disorder. Note that we have averaged here over ten configurations, an average over more 
configurations gave no better data. These two functions are related to each other by 

r=l 

which is exact for periodic boundary conditions. 

Due to numerical limitations the largest possible system size was L — 128. Here we found that the depinning 
transition takes place at he — 0.0029. We have analyzed the morphology of the domain wall for a driving field 
slightly above this critical field h — 0.003. Fig. (la) shows the height correlation function as a function of distance 
r for different times on logarithmic scales. For small r a linear behavior is observed, i.e. C(r, t) shows power law 
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behavior, while for large r saturation sets in at C-values which depend on time. For this behavior which is well known 
from kinetic roughening phenomena usually a dynamical scaling form 



C( r , r )=£( r )^_j (9) 

works well. Here £(r) oc r 1 / 2 denotes a time dependent correlation length and z = a/ 13 a dynamic exponent. The 
scaling function g{x) has the properties g(x) « const for x 1 and g(x) oc x 2a for x <C 1. The above scaling form 
contains two important limiting cases. The first one, 

C(r>e(r),r) oc r 2/3 , (10) 

describes the growing of the spatially uncorrclatcd interface fluctuations with time whereas the second one describes 
the growing of the spatial correlations: 

C(r<£(r),T) ocr 2Q . (11) 

Analyzing the height correlation function in the above limiting cases we get a — 0.68±0.01 for the roughness exponent 
and for the small time exponent (3 = 0.36 ± 0.02 (see inset of Fig. (la)). Note that we have performed simulations for 
different strength p of the random fields and found no dependence of the exponents on p. Due to the occurrence of an 
intrinsic width |TT| of the interface, which manifests it selves in a turning point of the height correlation function for 
small r, the scaling form Eq. (0) holds only for r > 3. The average step height of the wall, characterized by C(l, r), 
shows a small time dependence, see Fig. (la). If we fit C(1,t) to a power law, C(l,r) oc t 2k , as we have done in 
d = 1 + 1 [ p"l| we find a very small value for this exponent, k « 0.014, on a time interval which is smaller than the 
interval at which the scaling behavior Eq. ( |l0| ) is valid. Because of the smallness of this exponent and due to the 
small time interval where an algebraic behavior can be observed we think that it makes not much sense to extend 
the scaling law Eq. (|9|) as had to be done in d = 1 + 1 where this exponent is much larger |11|. Fig. (lb) shows a 
satisfactory scaling plot. 

Since in the limit u — ► oo the present soft-spin model goes over into the random-field Ising model (RFIM) jl3) it is 
tempting to compare with this model. The theoretical prediction of Grinstein and Ma |I3] for the RFIM is a — (5— d) /3 
while Ji and Robbins jl5| obtained from numerical studies of the RFIM in d = 3 the value a = 0.67 ± 0.03. Amaral 
et.al. H found that for the RFIM the prefactor of the KPZ-nonlinearity is zero or goes to zero at the depinning 
transition so that the RFIM should belong to the EW-universality class. The value of the roughness exponent we 
obtain is close to the above cited values. Note that Leschhorn [jl6| found for a solid-on-solid (SOS)-type modification 
Eq. (|l|) ct = 0.75 ± 0.02 significantly larger than our value. 

With the values for the roughness exponent a and the small time exponent (3 we obtain for the dynamic exponent 
z = a//3 — 1.89 ± 0.13. We consider this value as a confirmation of our dynamic scaling analysis for Eq. ([!]) presented 
earlier |D|] where we found z = 2 for all d and a = (5 — d)/3. Note that z » 2 was found previously also in d = 1 + 1 
for the model discussed here [ p"l| and again in d = 1 + 1 from a direct integration of the equation of motion Eq. (|l|) 

B 

The result for the dynamic exponent and therefore for the small time exponent f3 = a/z disagrees with the numerical 
results of Leschhorn for his SOS-model who found (3 = 0.475 ± 0.015 and z = 1.56 ± 0.06 in d = 3. The reason for 
the discrepancy of the numerical values of Leschhorn and others is not clear to us especially since Amaral et.al. H 
could show that Leschhorns SOS model can also be described by Eq. (|l|). A renormalization group study of Eq. ([j]) 
has been done by Nattermann et al. [jl9| who found z = 2 — 2e/9 in d = 5 — e. If one extrapolates this down to d = 3 
which, however, might be well outside the range of the linear e-expansion one obtains rather good agreement with the 
value of Leschhorn but not with our result. Still further work has to be done to clarify this situation. 
The situation described so far is only valid for systems near the depinning transition. Far away from it, i.e. for 
driving fields h 3> he the domain wall moves with a large velocity and the quenched random fields act as an effective 
white noise. Additionally, a nonlinear KPZ term is important H. Simulations in this regime are extremely difficult 
if not impossible since one has to go to very large system sizes. For small systems EW-behavior is observed and the 
crossover to KPZ-behavior occurs only for very large system sizes in d = 1 + 1 and for system sizes practically outside 
the range for simulations in d = 2 + 1 |2Q,^,[f8[ . 
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FIG. 1. a) The height correlation function C(r, r) for h = 0.003 at various times as indicated in Fig. (lb). The line 
represents the fit to C(r, r — > oo) oc r 2c * (Eq. (|l l[) ) . The inset shows the behavior of C(L/2, r)._The line represents the fit 



C(L/2, 



(Eq. 



driving field as in (a). 



b) Scaling plot of the height correlation function according to Eq. for r > 3 with the same 
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